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1. Introduction 

Rare decays of K-mesons as well as K^-K^^ mixing continue to play an important role in fixing 
parameters of the standard model (SM) and in constraining models of new physics. In the future 
Ek, the parameter describing indirect CP violation in kaon mixing, and the decays — > tt+vv, 
Kl ^ n'^vv will provide a decisive test of the SM and its extensions: they are highly sensitive to 
new physics [1] and their theory prediction is under good control for Ek and remarkably clean for 
K nvv. In the SM these modes are dominated by internal top quark contributions proportional 
to powers of Vf*Vtd and as such are suppressed with respect to generic new physics scenarios by 
the near diagonality of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. Also these modes can be 
calculated precisely using an effective theory framework. The matrix elements are extracted from 
Kit, decays for K — > 71 vv [2], and from the lattice for Ek [3]. This leads to an exceptionally clean 
prediction for the rare K —>■ nvv decays, while the recent and expected progress for the lattice 
calculations of Bk, the bag parameter for Ek, is quite remarkable. 



2. Structure ofK^ nvv at NLO 



The theoretical cleanness of the K — > nvv decays in the SM is related to the quadratic Glashow- 
Iliopoulos-Maiani (GIM) mechanism. Using A,- = and Xi = mj/M^ we can write the ampli- 
tude of the Z-penguin and electroweak box diagrams (Fig. 1) as 



X, {F{xt)-F{xu))+Xc{F{xc)-F{xu)) 



A4Vn^ + ^U^]ln"' 



Ml, J Mw \ r 'Mw 



(2.1) 



Here XtF{xt) is the top quai^k contribution, which is suppressed by five powers of the Cabibbo angle 
X = I Vus|> and XcF{xc) is the charm quark contribution. The contribution of soft internal up quarks 
is suppressed by Aq(^]-,/M^. 

Related to the quadratic GIM mechanism is the fact that the low-energy effective Hamiltonian 

^eff=^ ^ I (XcX\x,)+XrX{x,)){sLy^dL){Vli^y^Vli) (2.2) 



^/l 27rsin^% , 

involves only one single operator Qy = {slJ^ di){vLyp.VL) ■ The hadronic matrix element of the 
low-energy effective Hamiltonian can be extracted from the well-measured Kit, decays, including 
isospin breaking and long-distance QED radiative corrections. 

The anomalous dimension of the operator in (2.2) vanishes. Hence the top quark contribution 
in (2. 1) has no large logarithm and is calculated in fixed-order perturbation theory. The matching of 
the electroweak box and Z-penguin diagrams in Fig. 1 with internal charm quarks gives the charm 
quark contribution to the Wilson coefficient of Qy. Next, one matches Green's functions with 
internal W- and Z-Bosons and dimension-six current-current operators. The bilocal mixing into 
the dimension-eight operator rn^Qy - see Fig. 1 - resums the large logarithm in (2.1). The GIM 
mechanism cancels all loop contributions which do not carry an explicit charm mass dependence. 
Only when integrating out the charm quark higher-dimensional operators appear. The matching 
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Figure 1: Left column: The Z-penguin and electroweak box give the matching contribution to the charm and 
top sector. Middle column: The mixing of the current-current operators into Qy and the self-mixing of 
the current-current operators. Right column: Integrating out the charm quark produces m^Qv and subleading 
higher dimensional operators. 



onto 2v in Fig- 1 gives the dominant contribution to the branching ratio, while the contribution 
of the higher-dimensional operators can be computed together with the matrix element containing 
soft up quarks with the help of chiral perturbation theory ix^) [4]- 

After extracting the matrix element of 2v from Ki^ decays and summation over the three 
neutrino flavours the resulting branching ratio for K —>■ nvv can be written as [4, 5, 6] 

2 



TT+VV) 



k:+(1+Aem 
/ImAf 



2 



(2.3) 



Here k:+ = 0.5173(25) x 10-"'(A/0.225)^ and JCt = 2.231(13) x 10-i''(A/0.225)*^ contain higher- 
order electroweak corrections for the normalisation to the Kit, decays, and Aem — —0.3% denotes 
long distance QED corrections [6]. The top quarkX(;c,) = 1.464ib0.041, computed at two-loop in 
Ref. [7], gives the only contribution to ■S^{Kl K^vv), while its contribution to SS{K^ tt+vv) 
is 63%. The perturbative part of the charm quark contribution at NLO [7] is 

dNLO 



0.364 ± 0.036theory ± 0.009„,^ ± 0.009«^ 



(2.4) 



for A = 0.2255. The parametric uncertainty - see Ref. [8] for input parameters - is small compared 
to the theoretical error, which results from higher order corrections. In a ;^PT calculation [4] the 
contribution of higher dimensional operators and soft up quarks has been calculated to 5Pc,h = 
0.04 lb 0.02. Using Eq. (2.3), Eq. (2.4), and the input parameters of Ref. [8] results in: 

-11 

(2.5) 



(8.5±0.5sd±0.2id±0.6param)xlO 



(2.7±0.1sd±0.04id±0.4 



param. 



X 10 



11 
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Figure 2: The A5 = 2 Box-type diagram with internal up, charm, and top contributions is expressed as a 
sum of Box-type diagrams proportional to A/, X^, and XiXc respectively using the GIM mechanism. 



The subscript "sd", "Id", and "param" labels the perturbative, long-distance, and parametric uncer- 
tainties respectively. The parametric error is dominated by the CKM parameters, while 5Pc_„ gives 
the largest contribution to the long-distance uncertainty of the charged decay mode. 



3. Structure of £k at NLO 

The parameter 



measures CP violation in ^o-^o mixing via the ratio of the respective decay amplitudes of a Kl 
and a Ks decaying into a two pion state of isospin zero. Its generic structure shares features with 
the previously discussed K nvv decays: only one operator Qs2 = {sy^dL){sy^di) contributes 
dominantly below the charm quark mass scale and long-distance effects are power suppressed, this 
time because of CP violation. 

For the theoretical prediction it is useful to express Ek in terms of {K^\JPj^j"''^\EP) = lmKMl2, 
the matrix element of the A5 = 2 effective Hamiltonian, and write: 



Here the phase of Ek is 0e = 43.5(7)° [9] and t, = ImAo/ReAo ~ is the imaginary part divided by 
the real part of the isospin zero amplitude Aq = {Ks — > (7r7r)/=o)- The ratio Kg = \e^ / EK{<^e = 
45°, i§ = 0)1 encompasses the change of \Ek\ if the values 0e = 45° and = are used in (3.2), as 
has been done in most of the older analyses, instead of the exact values. The authors of Ref. [10] 
used e' /e to extract the value of JCg = 0.92 ± 0.02 in the SM. 
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a) b) 



Figure 3: Dimension 6 current-current and penguin operators mix at LO into Qs2 with a CKM factor 
proportional to ?^Xc in a). Integrating out the charm quark results in similar diagrams for the LO and NLO 
matching of the contribution proportional to XtXc and respectively. A sample diagram which is relevant 
to the LO evolution of gs2 is shown in b). 

The box diagram of Fig. 2a gives the leading contribution to the AS = 2 effective Hamiltonian 
and the parameter Mu- It is proportional to a sum of loop functions times CKM factors: 

9 2 
fU TTl 

£ A„,A„/(-^, = X^S{xr) + X'-S{x,) +Xr?icS{xr,x,) . (3.3) 
Ui.iije{u,c.t} '^W '^w 

After the GIM mechanism has been used to eliminate A„ = —A, — Xc it comprises the top quark 
contribution - proportional to (Fig. 2b), the charm quark contribution - proportional to 
(Fig. 2c), and the charm-top quark contribution (Fig. 2d) - proportional to XfXc- The resulting 
loop functions S{xi,Xj) = S{xi,Xj) —S{xi,0) —S{0,Xj)+S{0,0) and S{xi) = S{xi,Xi) are suppressed 
by the smallness of the quark mass m,- if Xi = mj/My^, is significantly smaller than one. This, 
together with the severe Cabibbo suppression of the CP violating top quark contribution, lets all 
three contributions compete in size for Sk'. 

Im {X^S{xr) + X^S{x,) + W{x„x,)) ^&{X'') + if (a^^) In + ^ (a^^) . 

(3.4) 

The diagram of Fig. 2a induces a large logarithm (Inm^/M^y) only for the chaim-top quark contri- 
bution: the lai^ge logarithm from the up-quarks in Fig. 2b is power suppressed by A^^^/M^, while 
the GIM mechanism cancels a potential Inmc/Miy between the diagrams with both one up and one 
charm quark and the diagram with only internal charm quarks. 

This can be reformulated in an effective theory language: the dimension-six penguin as well as 
the current-current operators, which have tree-level Wilson coefficients, mix only into the charm- 
top quark contribution, via the bilocal mixing in Fig. 3a, yet do not induce large logarithms times 
tree-level Wilson coefficients proportional to X} and A^^. QCD corrections do not change this 
picture but only induce the well known RGE effects for the AS = 1 effective Hamiltonian [11] and 
for the AS = 2 Operator Qs2 (Fig. 3b). A LO analysis of the charm quark and top quark contribution 
to Ek then requires a one-loop calculation both for the matching at /iiy, for the running, and for the 
matching of the charm quark contribution also for the matching at jic (Fig- 3a). This is contrary to 
the chaim-top quai^k contribution where a tree-level matching at pLw and pLc is sufficient at LO. 

After integrating out the charm quark the AS = 2 effective Hamiltonian reads 

=^eff = -AmI [X^^,,S{x,) + X'-^,,S{x,) +X^^aS{xc,Xt)\ b{^)Qs2 + h.c. + ... (3.5) 
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wiiere the QCD and logarithmic corrections are known at NLO and parametrised by T]cc = 1 .43(23) 
[12], T7ct = 0.47(4) [12], and % = 0.5765(65) [13]. The parameter ^(/i) is factored out such that 

3 , ,(^"|G52|^) 

BK=-b{^ y j ' (3.6) 

is a renormalisation group invariant quantity, which can be calculated on the lattice - see e.g. [14] 
andRef.[3] of this conference. Using S^f = 0.720 ±0.013 ±0.037 one finds for at NLO [10]: 

gNLO^^^ 7g^Q_25) , (3.7) 

where T7tt, T]ct, and T]cc contribute with 75%, 37%, and —12% respectively to the total value of 
£k, while 60% of the uncertainty is of parametric origin and 40% is of theoretical origin. The 
parametric error is dominated by the uncertainty in the CKM parameters, while the perturbative 
and non-perturbative uncertainties are comparable in size for the theory uncertainty. 

Finally note that J^^^^ also contains higher dimensional Operators in and current- 

current operators with up-quarks in ^g^J^^ as indicated by the ellipses in Eq. (3.5). At leading 
order in the l/Nc expansion only one higher-dimensional operator is present and its matrix element 
is estimated in [15] to result in a 0.5% enhancement of Ek- 



4. NNLO QCD corrections 

The NNLO calculation of Ek and K^vv aims at resumming all (^(a"ln""^(ju|,//i^^)) 

logarithms for Pc and T]ct and all {a^ In"^^ {|J.^, / jj,^)) for X{xt), T]tt, and rjcc- The theory uncer- 
tainty of Pc and rjct dominates the perturbative error for n^vv and Ek respectively at NLO, 
while the large theory uncertainty in T]cc is somewhat suppressed by the smallness of the charm 
contribution to Ek- 

A NNLO analysis for P^ and T]ct will reduce the theory uncertainties and comprises (/) the 
^{ocj) matching corrections to the relevant Wilson coefficients arising at Hw, (U) the &{(X^) 
anomalous dimensions describing the mixing of the dimension-six and the Qy and 2s2 opera- 
tors, (///) the i&{ccj) threshold corrections to the Wilson coefficients originating at Hb, and (iv) the 
^{ccj) matrix elements of the operators emerging at Hc- To determine the contributions of type (/), 
(///), and (/v) one must calculate two-loop Green functions in the full SM and in effective theories 
with five or four flavours. Sample diagrams for steps (/) and (/v) are shown in the left and right 
columns of Fig. 4. The contributions (//) are found by calculating three-loop Green functions with 
operator insertions. Sample diagrams with a double insertion of dimension-six operators are shown 
in the centre column of Fig. 4. The corresponding three-loop amplitudes are evaluated using the 
method that has been described in [1 1, 16]. A comprehensive discussion of the technical details of 
the matching, the renormalisation of the effective theory and the actual calculation is given in [5] 
for the calculation of Pc and will be given in [17] for the calculation of T]ct. The same techniques 
have also been used to reduce the uncertainties in the short-distance contribution to Kl — > H^H^ 
[18]. 

The aforementioned QCD calculation results for the input parameters of Ref. [8] in the fol- 
lowing value for Pc at NNLO: 

pNNLO ^ Q ^ 0.009theory ± 0.009^^, ± 0.009a, . (4. 1) 
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Figure 4: Examples of Feynman diagrams arising in the full SM (left column), describing the mixing of 
operators (centre column) and the matrix elements (right column) for the Z-penguin sector of TT+VV 
(upper row) and for Ek (lower row). Only the divergent pieces of the diagrams displayed in the centre column 
have to be computed, while the Feynman graphs shown on the left- and right-hand side are needed including 
their finite parts. 

Comparing these numbers with Eq. (2.4) we observe that the NNLO calculation reduces the the- 
oretical uncertainty by a factor of 4. Because of this and the improvement in the calculation of 
long-distance effects [2], unknown electroweak corrections could potentially dominate the theory 
uncertainty of the rare K —>■ Tivv decays. Even though a similar reduction of the error for Tjct is 
expected at NNLO in QCD [17] no electroweak corrections are needed for the present theoretical 
status of Sk- 

5. Electroweak corrections 

The NLO calculation of electroweak corrections for rare K nvv decays resums all LO and 
NLO logarithmic QED corrections and fixes the scheme electroweak input parameters, like sin^ dw, 
by an electroweak matching calculation. The function depends on the charm quark MS mass 
through the parameter Xc, conventionally defined as Xc = m^./M^. The point of fixing the input 
parameters can be exemplified by noting that the charm quark contribution is mediated by a double 
insertion of two dimension-six operators. This results in a contribution of i^{Gp) - the second 
power of Gf resides in x^ - plus electroweak corrections. Yet the leading result of Eq. (2.2) can 
only approximate the electroweak corrections for a specific choice of the renormalisation scheme 
for the pref actor of the charm quark contribution, expressed as a/ sin^ dw- While it is expected that 
using MS parameters renormalised at the electroweak scale would approximate the electroweak 
corrections best [19], only an explicit calculation can provide a definite result. We normalise all 
dimension-six operators to Op and replace the parameter x^ with the definition 

Xc = V2^^^^Gfm2(jU,) , (5.1) 
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Figure 5: Left column: Example of a diagram describing the NLO matching to a dimension-six opera- 
tor involving charm quarks and neutrinos (top) and of a diagram contributing to the NLO mixing of two 
dimension-six operators into Qy. Right column: Pc{X) as a function of at NNLO QCD (dashed dotted 
line), including LO QED (dotted line), and NLO electroweak corrections (solid line). The dashed line shows 
PciX) at NNLO QCD where the definition Xc = /M^, is used. 

which only at tree level equals the ratio mKnc) /M^. 

The NLO analysis of electroweak effects of Ref. [8] involves the calculation of one-loop 
matching corrections for the dimensions-six operators (top left diagram of Fig. 5) and QED cor- 
rections to the LO QCD operator mixing (bottom left diagram of Fig. 5), and the inclusion of QED 
effects in the expansion of the matrix elements at pL^. Note that the LO QED corrections start at 
ff{a^\vL'{jX^/}x]:)) while the first NLO electroweak correction is ^(a^ln(/i|,//i^)). This explains 
why Pc{x), which is plotted on the right column of Fig. 5 as a function of the parameter jic, receives 
corrections of similar size. Also the cancellation of the scheme dependence between the LO QED 
and the NLO electroweak contribution is clearly visible and we see that including the full elec- 
troweak corrections, Pc{^) is mildly increased as compared to the pure NNLO QCD. The number 
for the branching ratio then reads [8]: 

^ {K+ ^ K+vv) = (8.51+!^:^^^^^ ±0.20,„^,„,„a, ±0.36theory) X IQ-" . (5.2) 

The CKM parameters dominate the parametric uncertainty. The main contributions to the theory 
error stem from the uncertainty in 5Pc,u and Xt, where we used an error of 2%. In detail, the 
contributions to the theory error are (K-+ : 6%, Xt : 38%, Pc : 17%, 5Pc,u ■ 39%), respectively. All 
errors have been added in quadrature. 

6. Conclusions 

Rare K nvv decays and the CP violating parameter Sk are extremely sensitive to flavour 
violating new physics. The good control of long-distance contribution to these observables makes 
the calculation of NNLO QCD and sometimes even NLO electroweak corrections mandatory. Re- 
sults for NNLO QCD and NLO electroweak corrections for the charm quark contributions to rare 
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K decays have been published in Ref. [5, 8], while the NNLO calculation of the charm-top quark 
contribution is finished [17]. This, together with cuiTcnt [20] and future [21] progress from the 
experimental side, will increase the new physics reach of these observables further. 
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